Coherent structures in localised and global pipe turbulence 
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The recent discovery of unstable travelling waves (TWs) in pipe flow has been hailed as a signifi- 
cant breakthrough with the hope that they populate the turbulent attractor. We confirm the exis- 
tence of coherent states with internal fast and slow streaks commensurate in both structure and en- 
ergy with known TWs using numerical simulations in a long pipe. These only occur, however, within 
less energetic regions of (localized) 'puff' turbulence at low Reynolds numbers (Re = 2000 — 2400), 
and not at all in (homogeneous) 'slug' turbulence at Re = 2800. This strongly suggests that all cur- 
rently known TWs sit in an intermediate region of phase space between the laminar and turbulent 
states rather than being embedded within the turbulent attractor itself. New coherent fast streak 
states with strongly decelerated cores appear to populate the turbulent attractor instead. 



The transition to turbulence in wall-bounded shear 
flows is a classical problem that has challenged physi- 
cists for over a century. Some flows, such as that be- 
tween differentially-heated parallel plates or between ro- 
tating concentric cylinders, exhibit a smooth progression 
to increasingly complicated flows via an initial linear in- 
stability. Plane Couette flow and pipe flow, however, 
abruptly adopt a turbulent state. The problem is further 
complicated by changes in the spatio-temporal character 
of the observed flows at different flow rates. Pipe flow 
exhibits a quasi-stable localised turbulent 'puff' state as 
well as a globally turbulent 'slug' flow A finite am- 
plitude disturbance is required to trigger turbulence, the 
amplitude of which has been shown to depend critically 
on its 'shape' [2j]. Similarly, plane Couette flow exhibits 
localised spots of turbulence which may be either short 
lived transients or survive to arbitrarily long times Q. 
Plane Poiseuille flow, on the other hand, exhibits a linear 
instability but at flow rates well beyond those at which 
turbulence is typically observed. 

The discovery of exact travelling wave(TW) solutions 
in wall-bounded shear flows 0, 0, U [H, S G3 has spurred a 
flurry of excitement within the community and prompted 
the speculation that the states lie inside the turbulent at- 
tractor. Since TWs appear only to have a few unstable 
directions in phase space [1, EH, a turbulent flow trajec- 
tory is imagined to wander between these states, with 
the probability of the flow 'visiting' a particular TW de- 
termined by how unstable it is. The hope is that turbu- 
lent statistics are then predictable from an appropriately- 
weighted sum of all the relevant TW properties (see the 
reviews [HtDl)- 

Pipe flow has emerged as the main setting to con- 
firm this picture. Thus far, numerical results have been 
confined to unrealistically short periodic pipes G3 , 
and while experimental observations suggest evidence of 
TWs, important structural differences remain (compare 
figs 2E & F and 4 in [15] and fig 1 in [16]). In this Let- 
ter, we use direct numerical simulations in a sufficiently 
long pipe to capture real localised ('puff') and global 
('slug') turbulent states, in order to determine whether 
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FIG. 1: Skin friction coefficient C f = -(d z p) r e z t D/2pU 2 for 
a 50 D pipe where p is the fluid density. Triangles - laminar 
flow; Circles - puffs; solid Squares - spatially-inhomogeneous 
turbulence and open Squares - homogeneous slug turbulence 
(the transition is gradual). The lower straight line corre- 
sponds to laminar flow and the upper to the Blasius friction 
law which is initially overshot by the turbulent flow |4(. At 
Re = 2300 the flow alternates between a 1 puff state (Circle) 
and 2 or 3 puff states (Square with errorbars). 



TWs populate the turbulent attractor. The ability to ac- 
curately simulate spatially-inhomogeneous turbulence at 
transitional Reynolds numbers proves crucial in revealing 
the true position of the TWs in phase space. 

The non-dimensionalised governing Navier-Stokes 
equations for an incompressible fluid are 



dtu + u • + Vp • 



Re 



V-u = 0, (1) 



where the Reynolds number Re := UD jv (U is the mean 
axial flow speed, D is the pipe diameter and v is the 
fluid's kinematic viscosity). These are solved in a pipe of 
length 50 D across which periodic boundary conditions 
are imposed at constant Re (mass flux) in cylindrical 
coordinates (r, z) using the formulation and numerical 
resolution described in [17[. (At Re = 2800, the energy of 
spectral coefficients falls by at least 6 orders of magnitude 
with spectral order in either r, 6 or z.) 

Data was collated from simulations over a range of Re 
up to 3500 and over times of greater than 3000 D/U; 
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FIG. 2: Axial component of vorticity in (r, z)-plane, 25 D shown of 50 D computational domain; (a) slug turbulence at 
Re — 2800; (6) inhomogeneous turbulence at Re — 2400 (c) puff turbulence at Re — 2000. Cross-sections in (r, 9) show the 
axial flow relative to the laminar profile with fast streaks (light/white) and slow streaks (dark/red), contour lines each Q.2U; 
(d) m — 4 and m — 3 structures seen upstream and downstream of the trailing edge zte (flow is left to right); (e) sections from 
a puff where no clear structures are observed; (f) energetic section at Re = 2800, but resembling m — 5 structure; (g) exact 
solution with 3-fold rotational symmetry. 



the skin friction is a useful indicator of the flow response 
- see Fig. [TJ A localised 'puff' structure [l| of appar- 
ently stable length w 20 D is observed for Re < 2250. 
At 2250 < Re < 2500 the puff gradually expands while 
translating along the pipe, possibly dividing into multi- 
ple puffs, giving rise to an uneven patch of turbulence in 
which the turbulent intensity is spatially inhomogeneous 
(see Fig. [2]a-c). By Re w 2800, the turbulent intensity 
has become much more spatially homogenized indicating 
'slug'-like turbulence 

As a puff is spatially inhomogeneous, the search for 
coherent structures was conducted at fixed relative posi- 
tions, up- and down- stream, of the puff's steadily trans- 
lating trailing edge, zte(£)> which is itself characterized 
by a sharp jump in the streamwise velocity u z on the 
pipe axis [22[. The search focused upon the appearance 
of fast streaks near the pipe wall, thus correlations in the 
streamwise velocity were examined using the function 



C(0,z') 
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where ( • ) s indicates averaging over the subscripted vari- 
able, and u' z is the deviation from the time-averaged pro- 
file calculated for each z' = z — zte position in the puff. 
The projection function C m {z') = 2 (C(0, z r ) cos(m0))o 
was used to extract the signature of structures of az- 
imuthal wave number m. Experience of examining flow 
structures indicated that a 'good' correlation is achieved 
for C m (z f ) larger than 0.1, as indicated by Fig. [3l which 



shows the correlation results for the puff snapshot of Fig. 
[2]c. The magnitude of the correlation measures are rel- 
atively large at the positions indicated in Fig. [3] despite 
not being located at zte, where the turbulent intensity 
is greatest (u' z is largest) for the puff. Cross-sections in 
(r, 6) of the flow field are shown in Fig. [2]rf with lines in- 
dicating their position. For comparison purposes, cross- 
sections for another puff snapshot are reproduced in Fig. 
[2]e where C m (z') is less than 0.1 for all z' . Particularly 
for the plots upstream of the trailing edge, the similar- 
ity to known TWs is remarkable where individual slow 
streaks are also reproduced in the interior (compare, for 
example, with Figs. 9a and 13(lower left) in [9j]). 

The probabilities of finding a correlation greater than 
0.1 at different parts of the puff are plotted in Fig. [H 
From around z' = —D to z' = +5D fast streak struc- 
tures of m = 3 and m = 4 are seen approximately 10- 
15% of the time which is in good agreement with the 



frequencies observed in [11|, [14|, [18[ . The appearance of 



coherent fast-streaks, however, does not necessarily im- 
ply an observation of a TW. Of the coherent fast streak 
structures found, those which look most like TWs (in 
terms of fast and slow streaks) are found away from the 
most energetic regions in puff turbulence. The distur- 
bance energy at the trailing edge itself (Fig. [5]) is far 
too high to be compatible with any known TW at the 
same Re with the roll energy, in particular, an order 
of magnitude too large. Cross-sections at zte also ex- 
hibit small-scale structure uncharacteristic of TWs (see 
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FIG. 3: Instantaneous correlations for the puff snapshot of 
Fig. [2]c at positions relative to the trailing edge z — z — 
zte • In the upper plot goes from to 2n vertically, contour 
intervals of 0.1. The marks below indicate the positions of 
the cross- sect ions of Fig. [2]g?. 



Fig. Wo)- Upstream and downstream, at approximately 
z' = —2D and z' = where we find well formed 

coherent structures, the magnitudes of streak and roll 
energies are both consistent with TWs which all have a 
characteristically small roll-to- streak energy ratio. 

At higher Re, the turbulence becomes delocalized as 
the puffs expand to invade the whole flow. The strongly 
turbulent region around the trailing edge also lengthens 
to swallow up the weaker relaminarization zones so that 
the flow spatially homogenizes. Correlations can then be 
further averaged over the pipe length C m := (C m (z)) z 
where u' z is now taken as the deviation from the full- 
space-and-time-averaged profile. Fig. [H shows a clear 
trend in which the preferred structures are gradually de- 
creasing in scale from m = 3 and 4 streaks at Re = 2000 
to m = 4 and 5 streak structures at Re = 2800. A 
typical correlation episode in split puffs at Re = 2400 
is shown in Fig. [6] using three snapshots AD/U apart. 
There are two coherent structures simultaneously present 
in the 25 D section shown which translate with phase 
speeds of ~ (1 ± 0.1 )U. As in the puff, however, tran- 
sient signatures of TWs are found to occur in the less 
energetic regions of puff turbulence. 

At larger Reynolds numbers, Re > 2800, the turbulent 
intensity in slug turbulence is uniformly high everywhere, 
being comparable to that at the trailing edge of a puff (in 
units of pv 2 at each cross-section), and is never as low 
as that of the TWs (see Fig. [5]). Although fast streak 
structures are still observed, they are of too high energy 
to be associated with known TWs. The cross-section 
shown in Fig.EJf (C5 > 0.1) from a slug is typical, where a 
large and strongly retarded central core dominates, with 
only narrow fingers extending towards the gaps between 
the intense fast streaks at the wall [23[. This retarded 



core feature of the coherent structures found at Re > 
2800 can now be appreciated as a significant problem 
in previous qualitative comparisons between slug cross- 
sections and known TWs [HI, [l6[ • It is more likely that 
these coherent structures instead point to the existence 
of other types of exact, more highly nonlinear solutions. 

Our results suggest that the known TWs populate an 
intermediate region of phase space between the laminar 
and fully turbulent phases, rather than the turbulent part 
of phase space itself. The dynamical importance of the 
TWs is therefore in the transition-to-turbulence process 
where the fate of an initial disturbance is of concern 
rather than in characterizing the established turbulent 
state. The puff provides a simple illustrative example 
of this picture of TWs sitting between the laminar and 
turbulent states. At Re = 2000, a puff travels at only 
« 90% of the bulk velocity so, on average, fluid passes 
through it. Far upstream (t — > —00), the fluid 'trajec- 
tory' starts at the origin (laminar state), passes through 
the TW-region of phase space as t increases, to reach the 
fully turbulent region near the trailing edge. On leaving 
here as time increases, it passes back through the TW 
region to the origin as it relaminarises far downstream 
(t — > +00). Consequently, TWs are only visited just up- 
and down- stream of the trailing edge: Fig [5] is a good 
2-D representation of this process (where z plays the role 
of t and a phase space norm based on the streak and roll 
energies is implied). 

Our results confirm emerging evidence that lower 
branch TWs lie strictly between the laminar and tur- 
bulent states in phase space H, 0, 0, Hlj]. However, 
the fact that upper branches of the known TWs also have 
too low energy to be associated with the turbulent part of 
phase space is a surprise. It is quite plausible that the co- 
herent structures observed so far for Re > 2800 and char- 
acterised by an outer ring of fast streaks together with 
a strongly decelerated core represent a more-energetic 
branch of TWs which is embedded in the turbulent at- 
tractor. The fact that their roll-to- streak energy ratio is 
so much larger than for currently known TWs suggests 
that they may have a different sustaining mechanism. 
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FIG. 4: Probability at different parts of a puff of a 'good' correlation Cm(z') > 0.1 (see Fig. [3} at Re ■ 
of a correlation C m > C at any given point in the flow for Re = 2400 (middle) and 2800 (right). 
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^ 2 7r^ m _, / (\u mr \ 2 + \u' m9 \ 2 )rdr in units pv 2 . Here u = 



j-tn^O J ^ 1 '" 1 1 \ U m I 

(u' r ,u' e ,u z ) is the deviation from the laminar profile, for eas- 
ier comparison with the TWs. Energies for 2-, 3- and 4- fold 
rotationally- symmetric TWs [12] are shown using lower thick, 
middle thin and upper thick red solid lines respectively (the 
closed loops are produced by the finite continuum of TW axial 
wavenumbers which can exist at Re — 2000). The pink short- 
dotted line 'cloud' in the top right hand corner corresponds 
to a slug at Re = 2800. 
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